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We compute the mass spectrum for strange/charmed baryons in the partially quenched approach 
using Nf — 2 twisted mass QCD configurations. We investigate two main issues: the size of 
lattice artefacts using three values of the lattice spacing (the smallest of which is approximately 
0.05 fm) and the dependence of baryon masses on meson (or quark) masses. We thus perform a 
global fit in order to extrapolate simultaneously to the continuum limit and to the physical point. 
We estimate the masses of Q.sss, S^/js, A^ii^, Slccc, ^dcc, \idc- 
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1. Introduction 

Simulations with two light degenerate sea quarks {Nf = 2) and including also the strange sea 
quark {Nf = 2 + 1) are nowadays standard. The ETM Collaboration has generated a substantial 
sample of Nf = 2 ensembles at four values of the lattice spacing (ranging from 0. 1 to 0.05 fm), 
several values of the light sea quark mass and several physical volumes. Using these ensembles 
one can study the cut-off effects on observables and the insight gained provides valuable input for 
the choice of parameters for the Nf = 2 + 1 + 1 (i.e. including both the strange and the charm sea 
quarks) simulations under production. Preliminary results using Nf = 2 + 1 + 1 simulations have 
been presented at this conference 

In the present study we therefore use Nf = 2 ensembles with a partially quenched setup in 
which the strange and charm quarks are added only as valence quarks. For heavy quarks the 
compton wavelength of the associated heavy-light meson is small compared to present attainable 
lattice spacings which means that cut-off effects can be large. The charm quark mass is at the upper 
boundary of the range of masses that can be simulated at present for the coarsest lattice spacing 
used in the continuum limit extrapolation {a ~ 0. Ifm for which m^a < 1). In order to safely control 
this extrapolation it is thus important to asses the size of lattice artefacts affecting the observables 
of interest. 

Our goal is to extend the study of Ref. by including a finer lattice spacing a ~ 0.051 fm. 
We would like, in addition, to compute the low-lying spectrum of strange and charmed baryons. 
In this contribution we present preliminary results for the masses of the strange baryons Q-sss, '^dss, 
Kiids and the corresponding charmed baryon obtained by substituting the strange quark with the 
charm quark {Q-ccc, "^dcc, \idc)- Preliminary results for the low-lying strange baryon spectrum with 
Nf = 2 + 1 + \ gauge configurations [^] and for the spectrum of static-light baryons with Nf = 2 
configurations [Q] have also been presented at this conference. 

2. Setup 

The lattice discretization used for the doublet of degenerate light quarks is Wilson twisted 
mass QCD at maximal twist whose action reads (in the twisted basis) 

^hT° = (y (Vm + V;;) - ^V;V, +M„ + /75 T3 M/) Z/W (2.1) 

where V^, V* are forward and backward covariant derivatives, Mcr is the Wilson critical mass and 
/i/ is the light quark mass. 

The strange and charm (which in the following are referred to as "heavy") quarks are added 
here only as valence quarks a la Osterwalder-Seiler and their action reads 

X h=s ^ ^ ^ 

where \is and \ic are the strange and charm (valence) quark masses. In order to remove the deter- 
minant of the strange and charm quarks, ghosts have to be added correspondingly. Concerning the 
gauge actions, ETMC uses the tree-level Symanzik improved gauge action. 

The main advantage of this regularization with respect to the standard Wilson one is that the 
spectrum and the matrix elements extracted from correlation functions are automatically Oia) im- 
proved [^. The drawback is that parity and isospin are explicitly broken by O(a^) lattice artefacts 
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and are recovered only in the continuum limit. Here we use ETMC configurations generated at 
three values of the lattice spacing a G {0.051,0.064,0.080} fm and physical volumes L ~ 2.0-^2.4 
fm (the scale has been set through fj^ in Ref. Light sea quark masses correspond to pion 
masses Mji € [290, 520] MeV while partially quenched valence strange/charm quarks correspond 
to heavy-hght meson masses Mk G [520,710] MeV and Md G [1.80,2.40] GeV. In all we have 40 
different combinations {MjiMhi)- In order to combine data at different lattice spacings we express 
the value of the masses in units of ro [||]. For the three lattice spacings considered here the values 
ro/a G {8.36,6.73,5.36} are taken from Ref. ^ 

3. Numerical results 

An important issue in our study is the dependence of the baryon masses upon the "heavy" 
quark mass /i/, in the strange and in the charm region. At a = 0.080 fm and Mj^ ~ 340 MeV this 
dependence is shown in Fig. |l|. 
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Figure 1: The dependence of the octet and decuplet baryon masses on /i/,. Dashed lines corespond to linear 
fits performed in the strange region. 

From Fig. [I] we observe that baryon masses depend linearly on jih both in the strange and in the 
charm region but with two different slopes. This behavior will be further discussed in what follows. 
For what concern meson masses, in the case of the Kaon we observe a dependence M\ oc in 
agreement with the fact that the Kaon can still be considered a pseudo Goldstone boson. For 
the D meson instead we observe a dependence Mu oc jXf^ as predicted by heavy quark effective 
theory (HQET), with no evidence of \/pih term. In the following we will consider the functional 
dependence of baryon masses upon Mji and M/,/ because this allows to extrapolate to the physical 
point without knowing the values of the renormalized quark masses. The observations above imply 
that baryon masses depend quadratically on Mk in the strange region while depend linearly on 
in the charm region. 

From Fig. [I] and Fig. ^ it is also evident that the splitting between / = 1 /2 and J = 3/2 states 
(£/£* and S/S*) clearly diminishes with the increase of /i/,. In quark models, this observation 
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Figure 2: E/E* and E/E* splittings as function of /X/,. 
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is explained thanks to the fact that the spin-spin coupling part of the q — q potential is inversely 
proportional to the masses of the two quarks In HQET, the splitting of baryons containing one 
heavy quark (e.g. the ILudh/^ldt) proportional to 1 //!/,. 

Hadron masses Mu are extracted from the two point correlators Ch(0 = ^x(^(^'^)^^(0,0)) 
of the corresponding interpolating operators H at large time distances. The interpolating operators 
H are those of Ref. and to improve their overlap with the ground state we apply Gaussian 
smearing and use APE smearing for the links that enter the hopping function. At large Euclidean 
time separation the value of the hadron mass can be extracted by fitting the effective mass defined 



hyM^{t) = \\n 



Cnjt) 
a C}i{t+a) 



to a constant. 



It turns out that the statistical error on M^{t) for the strange baryons grows faster in time than 
in the case of the charmed baryons. In the case of the Q.sss and Q.ccc this is illustrated in Fig. ^. It is 
easy to show that the statistical error on M^^^^ (t) is 
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Figure 3: Left: M^^Jr). Right: M^^Jt). 
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Figure 4: Mn as function of M/,/. The dashed Hne is an interpolating form between the s and c region. 

where M/-/, is the mass of the hh meson made of an heavy and an anti-heavy quark. This phe- 
nomenon is then probably explained by the fact that the gap Ajj^^^ = Mq^^^ — \Mcc has a smaller 



value than A^j^,^ = Mq^^^ — jM^s- At the physical point M^^,^ 
r]s meson would have a mass Mj., ? 



1672 MeV while the unphysical 

ss = rjs meson wouia have a mass M^s ~ y 2M|- - sa 690 MeV [§] and the gap Aq^^^ ^ 640 MeV. 
In the charm case instead, the preliminary prediction from the present work gives Mq^^^ ^ 4730 
MeV while the cc meson can be identified with the rjc meson which has a mass Mf,^ = 2980 MeV. 
The gap Aq^^^ 260 MeV is therefore sensibly smaller than in the strange case. Presumably, this 
fact remains true for the values of the meson masses we have in our simulations but we still need 
to check numerically this conjecture. 

3a X dlld ^^(^(^(^ 

In Fig. ^ we present all the 40 data points for the Q mass fitted to a functional form which 
interpolates between the strange and the charm region. This plot already shows the smallness of 
lattice artefacts. The functional form reduces, in the strange region, to the form Mq = Mq +AM]j. + 
BMIj. In the charm region it reduces instead to Mci = D + EM^ + FMhi. The two forms fit well 
the data: using 13 data points in the strange region we obtain Xdot ~ 1-56; using 27 data points 
in the charm region we find Xdof ~ Lattice artefacts are visible in the strange region and 
the inclusion of a term Aqo^ to the functional form above lower the x^o f fro^n 1.56 to 0.92. For 
charmed baryons we do not see any cut-off effect. 

This is due however to the choice of studying the behaviour of baryon masses as function of 
meson masses. Had we chosen to study their dependence upon the renormalized quark masses /i/ 
and }Xh (obtained from the bare masses by multiplying them by taken from Ref. [|l0|]) we would 
have immediately remarked the presence of lattice artefacts, at least in the charm region. In this 
region, a fit to the form =D + E^i + F^h (i e- not including lattice artefacts) is not sufficient 
and gives a huge X^^f- In order to obtain a reasonable Zdof ~ ^-^^ "'^^ needs to add lattice 
artefacts (both /X/, -independent and /i/j -dependent). Fig. || shows the data points together with the 
curve obtained by plotting the fitting function after setting a = 0. 

Lattice artefacts are instead hardly visible in the strange region. Here, a fit to the form 
Mqi = A + B}Xi+B}Xh gives a reasonable Xdof = ^-^^ ^^'^ adding lattice artefacts (dependent or 
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Figure 5: Mq as function of /X/, in the charm region. The dashed line is the continuum Hmit obtained by 
plotting the fitting function described in the text after extrapolating to a = 0. 

independent on the quark masses) does not improve the fit. 

We remark that Mq, depends very mildly on Mj^ and therefore the extrapolation to the physical 
Mjc seems not to pose any problem. By interpolating also to the physical value of Mk we get the 
result Mn,„ = 1.86(20) GeV which is consistent with the analysis in Ref. [Q] but still 10% larger 
than the experimental value. Due to the previous considerations and the analysis performed, this 
discrepancy seems not to be related to the continuum limit extrapolation or to the extrapolation in 
the light quark mass. Extrapolation to the physical {M„,Md) point gives the prediction Mjj,^^ = 
4.73(40) GeV (the experimental value is not known) and the ratio Mq^^JMq^^^ = 0.393(54) 

3.2 Auch and A^dc 

In the case of the A baryon, the dependence on Mj^ is much stronger than in the previous case 
and the inclusion of the term proportional to is crucial and reduces the xl^^ of a factor ~ 0.5 in 
both the strange and the charm region. Lattice artefacts are hardly visible and the functional forms 
we have used to fit are Ma = Mq + AM^ + BMli + CM^ (in the strange region) and Ma = D + EM^^ + 
FMhi + GM^ (in the charm region). Of course the inclusion of chiral logarithms would affect the 
extrapolation to the physical point. For these preliminary results we have however neglected them 
and performed only a rough fit using the forms written above. 

By extrapolating to the physical {Mj^,Mk) point we obtain Ma„^, = 1.20(10) GeV which has 
to be compared with the experimental value M^^"^^ = 1.116 GeV. By extrapolating to the physical 
{Mji,Mj)) point we have Ma„^^ = 2.24(18) GeV which is in good agreement with the experimental 
value M'^P =2.286 GeV. 

3.3 Zeiss and S^cc 

Twisted mass QCD breaks explicitly isospin symmetry and thus S°„. and ST (or equivalently 



S++ and S^^.^) are not degenerate. We thus preform a combined fit of both S j^^ and S^^.^ data with 
the form M^io.-} = Mq +AM1 + BM^, + CM]. + A{o where the coefficients A{o 
for the two sets of data. Analogously we perform a combined fit of both E 



-} 

■:r+ 



are different 



and EJ^.^ data with 
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the form M;^{++,+} =D + EM^ + FMhi + GM| + D^_^__^_ _^_'^a^. The dependence on M„ of the mass of 
the doubly charmed S turns out to be considerably less pronounced than for the standard (strange) 
E. The coefficient A_ is substantially smaller than Aq and in the charm case is compatible 
with zero and can be removed from the fit function. Fits work well and in the continuum limit, 
at the physical point, we get M^j,, = 1.37(12) GeV (to be compared with M't'^J' = 1.32 GeV) and 
Mej^^ = 3.52(25) GeV (in perfect agreement with m£P = 3.52 GeV). 

4. Conclusions 

In this preliminary study we have shown that, when baryon masses are analyzed as function 
of meson masses, lattice artefacts are always small and in some cases (notably in the charm region) 
hardly visible. They are instead clearly visible when baryon masses are analyzed as function of 
quark masses. As expected lattice artefacts are larger in the charm region, where they increase 
proportionally to Hh- The chiral extrapolation in the light quarks confirms to be critical and a 
term of order M| is needed for both S and A (it is particularly evident in this last case). Mq^^^ 
is still 10% larger than the experimental value and the source of this discrepancy seems not to be 
related to the continuum limit extrapolation or to the extrapolation in the light quark mass. Further 
investigations are needed to clarify this issue. Results for Ms^,,, ^A„d,, ^s^^c ^^'^ ^Kdc nicely agree 
with the experimental values. We have moreover obtained a prediction for Mq,^^^ = 4.13 (40) GeV. 
We are computing all the correlation functions needed to extract the whole low-lying spectrum 
of strange/charmed baryons. A complete analysis, including a more careful assessment of both 
statistical and systematic errors, will be performed in the near future. 
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